Applicability of Modified Effective- Range Theory to positron-atom and 

positron-molecule scattering 
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We analyze low-energy scattering of positrons on Ar atoms and N2 molecules using Modified 
Effective- Range Theory (MERT) developped by O'Malley, Spruch and Rosenberg [Journal of Math. 
Phys. 2, 49f (1961)]. We use formulation of MERT based on exact solutions of Schrodinger equation 
with polarization potential rather than low-energy expansions of phase shifts into momentum series. 
We show that MERT describes well experimental data, provided that effective-range expansion is 
performed both for s- and p-wave scattering, which dominate in the considered regime of positron 
energies (0.4 - 2 eV). We estimate the values of the s-wave scattering lenght and the effective range 
for e+-Ar and e"'"-N2 collisions. 



Electron (and positron) scattering on atoms at very 
low energies is dominated by polarization forces. Mod- 
ified effective range theory (MERT) was developed by 
O'Malley, Rosenberg and Spruch [J, Q for low-energy 
scattering of charged particles on neutral polarizable sys- 
tems in general. O'Malley applied MERT for elec- 
tron scattering on noble gases, in particular in the re- 
2pn of Ramsauer- Townsend minimum, using early total 
and momentum transfer [Slexperimental cross sec- 
tions. Haddad and O'Malley [g used three parameter 
MERT fit for s-wave phase-shift in electron-argon scat- 
tering and Ferch, Granitza and Raith 7] - for Ramsauer 
minimum in methane. Higher-order terms in MERT, 
resulting from short-range components of polarizabitity, 
were introduced by Ali and Eraser Jql- MERT analy- 
sis for Ne, Ar, Kr up to 1 eV was carefully revisited by 
Buckman and Mitroy |2l who used five-parameter fit for 
s-wave and p-wave shifts. 

Applicability of MERT to low-energy positron scatter- 
ing was already hypothesized in |2]- However, first mea- 
surements of total cross sections for positron scattering 
at low energies on noble atoms come only from seventies 

caini. 

The most systematic data for noble atoms, extend- 
ing down to 0.3 eV were done in WSU Detroit lab, us- 
ing positrons from a short-lived C^^ radionuclid, with 
about 0.1 eV energy resolution [Il|. Those data indi- 
cated clearly a rise of the cross section in the zero en- 
ergy limit in gases, like He, Ar, H2, Kr, X e, CO 9 , see 
|ll|. Unfortunately, subsequent experiments IZ, 13] used 
Ne^^ source and thick W-vanes positron moderator, thus 
worsening energy resolution and not allowing to make re- 
liable measurements below 1 eV. To gain in signal, large 
apertures and strong guiding magnetic fields were used, 
leading to underestimation of cross sections - some data 
showed even a fall in the limit of zero energy for highly 
polarizable targets, like CgHg [lj|- 

Only two of the most recent set-ups reached energies 
below 1 eV with good signal-to-noise ratio. In San Diego 



annihilation rates in Ar and Xe were measured, showing 
a steep rise below 1 eV [l^ . In Trento total cross sections 
in Ar and N2 were measured [1^ with angular resolution 
better by a factor of 30 than in some previous experi- 
ments jl2|]. Both laboratories confirm the early observa- 
tions from WSU Detroit on the rise of positron cross sec- 
tions in the zero-energy limit. Such a rise is also predicted 
by ab-initio theories [13, see [la| for detailed compari- 
son. A phenomenological attempt to apply MERT-like fit 
for low-energy cross sections in benzene and cyclohexane 
was done by Karwasz, Pliszka and Zecca IS]. 

In the present paper we apply MERT to positron total 
cross sections on argon and nitrogen, using recent ex- 
perimental data from Trento To|. We use MERT model 
based on direct solution of Schrodinger equation with po- 
larization potential as originally proposed by O'Malley, 
Spruch and Rosenberg jl|. Differently from earlier works, 
for the p-wave phaseshift we consider not only the po- 
larization potential but contribution from a general-type 
short range interaction. This allows us to extend the 
MERT applicability for positrons to energies above 1 eV. 
A clear indication on importance of p-wave scattering in 
this energy range comes from recent differential cross sec- 
tions measurements in argon JI3 • The present model in- 
troduces a second MERT parameter for the p-wave shift 
thus developing an approximation with two parameters 
for both s and p-wave phaseshifts. The first parameter is 
to be interpreted as a scattering length and the second as 
an effective range. We compare our MERT model with 
the expansion into momentum series valid at low energies 
and with ab-initio theories ]VA l20l| . 

Let us briefly review the effective-range expansion for 
1/r^ interaction. We divide the interaction potential be- 
tween charged particle and neutral atom into the long- 
range part: Vp(r) — — ae^/(2r^) with a denoting the 
atomic polarizability and e the charge, and the short- 
range part: Vs{r) describing forces acting at distances 
comparable to the size of atoms. In the relative coordi- 



nate the motion of particles is governed by 
r92 2d 1(1 + 1) {R*f 



Qj.2 



dr 



*W = o, (1) 



where "^{r) denotes the radial wave function for l-th par- 
tial wave, hk is the relative momentum of the particles, 
R* = ^J ae^p.jW- denotes a typical lenght scale related 



with the 



interaction, and /i denotes the reduced 



mass. In Eq. |^ we do not include V^(r\ which is re- 
placed by appropriate boundary condition at r ^ 0. The 
Schrodinger equation with polarization potential can be 
solved analytically jll l2lll2^ (see Appendix for details). 
At small distances (r ^ _R*) behavior of ^(r) is governed 

by 



*(r) 
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where is a parameter which is determined by the short- 
range part of the interaction potential. For r ^ i?*, $(r) 
takes the form of the scattered wave 



*(r)' 
with the phase shift rn 



1 
fcr 



sin(fcr — l^ + r]i), 



tan 77; 



m^ — tan S^ + B tan (5(m^ — 1) 
tan (5(1 - m2) 4- B(l - m^ tan^ 5) '' 



(3) 



(4) 



where we use similar notation as in Ref. [l|. Here 
B — tan(0 -I- l^), S = ^{v — / — |), and m and v are 
parameters obtained from the analytic solution of the 
Mathieu's differential equation (see Appendix). To in- 
troduce effective range we expand B around zero energy: 
B{k) = B{0) + ^RoR*k'^ + ... Jj. The second term can 
be intepreted as correction due to the finite range of the 
interaction, with Rq representing the effective range. 

In the zero-energy limit expansion of MERT in series 
of momentum k is useful. In the particular case of / = 
0, B{0) can be expressed in terms of s-wave scattering 
length Us'. B(0 ) = —R* /as, and expansion of cot?7o at 
fc = yields |lLl2^ 



qcotrio{q) 



3^'^ 



3a 



In 



1)'^ 



Rq 



2{R*y 



9a 3a2 9a^ 3a^^^' 



(5) 



where a = Ug/R*, q ~ kR* and ip{^) denotes the 
digamma function 24]. We apply similar procedure for 
p-wave. In this case, however, we expand directly tan 771 
given by Eq. Q 
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Here, 6 = B{Q) for Z = 1, and i?i denotes the effective 
range for p-wave. For higher partial waves we retain only 
the lowest order term in fc, which is sufficient to describe 
the scattering in the considered regime of energies 
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Let us turn now to positron scattering. We compare 
the total cross-section measured in experiment for Ar and 
N2 [13 with predictions of the theoretical model based 
on the effective-range expansion. In our approach the ef- 
fects of the short-range potential are included both for s- 
and p-wave giving the leading contribution to the scatter- 
ing in the considered regime of energies. Thus, our model 
contains four unknown parameters: the scattering length 
a and the effective range Rq for s-wave, and the zero- 
energy contribution 6 — B{0) and the effective range i?i 
for p-wave. For the investigated regime of positrons ener- 
gies, q = kR* can take values larger than unity, therefore 
for s- and p-wave we do not use expansions (01-® valid 
for (7 < 1, but we rather applied the initial formula © 
for the phase shift, performing only finite-range expan- 
sion for the parameter B. In this case, values of f and 
777, has to be evaluated numerically, using the approach 
described in the Appendix. 

For the calculations we use recent experimental values 
of the polarizability: a = 11.23ao^ and a = 11.54ao^ 
(atomic units), for Ar and N2, respectively [23. Table HI 
contains values of the characteristic distance R* and the 
characteristic energy E* — h^ / {2^R*'^) for the polariza- 
tion potential, and the values of four parameters: a, 6, 
-Ro: Ri which were determined by fitting our model to 
the experimental data. 

In the case of N2 the size of the molecule scaled by 
R* is much larger than for Ar, therefore we restricted 
our effective-range analysis to lower energies, fitting the 
model to experimental data with E < Q.%E*. In this 
regime the contribution of the effective-range correction 
in p-wave is rather small, and one does not get reliable re- 
sults for this parameter from the fitting procedure. Thus, 
for N2 we considered only three parameters a, 6 and i?o 
accounting for effects of the short-range part of the po- 
tential. 

Fig- n shows the experimental data for the total scat- 
tering cross-section for Ar as a function of positron col- 
lision energy. They are compared with: the MERT the- 
oretical curve which best fits the experimental data, its 
low-energy expansion given by Eqs. ©-(TI), and the re- 
sults of McEachran et al. [l3|- The total cross-section 
is presented in units of i?*, while the energy is scaled by 
E* . In the inset we additionally present contributions of 
the s and p waves to the total scattering cross-section. 
Similar results but for the scattering of positrons on N2 
are illustrated in Fig. 12 

We note the good agreement between our model and 
the results of McEachran et al. [I3 at low energies. 
Obtained value of the scattering lenght ag = -~5.5 8an 
agrees well with the calculations of McEachran et al. [13 



R*{ao) £*(eV) a^/R* 



Ro/R* Ri/R* 



100 



Ar 3.351 
Na 3.397 



1.211 
1.179 



-1.665 
-2.729 



-5.138 
-12.65 



0.3165 2.281 
0.8186 — 



TABLE I: Characteristic distance R* , characteristic energy 
E* and four fitting parameters: as (s-wave scattering length) , 
7?o (s-wave effective range), b (zero-energy contribution B(0) 
for p-wave) and Ri (p-wave effective range) for Ar and N2. 
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FIG. 1: Total cross-section for the scattering of positrons 
on argon versus the energy. Depicted are: experimental data 
(squares), the theoretical fit based on effective-range expan- 
sion (solid line), its low-energy part given by Eqs. lO-Qi 
and the theoretical results of McEachran et al. (dashed line). 
The inset shows in addition the s- and p-wave cross-sections. 
Data are scaled by the the characteristic distance R* and the 
characteristic energy E* of the polarization potential. 



(— 5.30ao), and Nakanishi and Schrader |23| (— 5.09ao). 
A somewhat worse agreement in N2 can partially result 
from poorer statistics of experimental data. Also ab- 
initio theoretical calculation in N2 show a big spread in 
determination of as, see for instance ^^. Both for Ar 
and N2 the model shows importance of p-wave scatter- 
ing above 0.7-0.8 eV. The s-wave effective range for Ar 
amounts to 1.06ao while in N2 to 2.78ao . We recall the 
" size" of the N2 molecule by the experimental determina- 
tion of the maximum HOMO density along the molecule 
axis which is about 2.3ao [23. Finally from Figs. Handle 
we observe that expansion into momentum series ©-0 
works only at very low energies below O.leV. 

In conclusion, we performed MERT analysis using an- 
alytical solution of Schrodingcr equation for polarization 
potential and we apply it to positron Ar (and N2) scat- 
tering up to 2eV. The scattering length in Ar agrees well 
with other predictions and the effective range (for s-wave) 
is 1.06ao. More experiments are needed at low energies 
to validate better the effective range parameters. 
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FIG. 2: The same as Fig.Qbut for the scattering of positrons 

on N2. 



APPENDIX A 

To solve radial Schrodinger equation 1^ we substitute 
r — \/R*e^^/\/fc and ^(r) = ip{r)^yR* /r, which yields 



the Mathieu's modified differential equation |24 12 



dz 



-j-j- [a-2qcosh2z]tlj^0. 



(Al) 



where a = {I + ^)'^ and q — kR* . Two linearly inde- 
pent solutions M(z] and T{z) can be expressed in the 
following form |22, |23| 



Mu{z) = Y. (-l)"c„(i^)J2„+.(2Vgcoshz),(A2) 

n— — 00 

00 

T.{z) = Y. (-l)"c„(i.)y2„+.(2V^coshz),(A3) 



which defines them for z > 0. Here, v denotes the char- 
acteristic exponent, and Jy{z) and Yu{z) are Bessel and 
Neumann functions respectively. Substituting the ansatz 
(|A2p and (|A3I) into (|A1|) one obtains the recurrence rela- 
tion: 



[(2n + vY - a] Cn + g(c„-i + c.„+i) = 0, 



(A4) 



which can be solved in terms of continued fractions. 
To this end we introduce /i+ = c„/c„_i and h~ = 
c-n/c-n+i for n > 0, which substituted into l|A4|l gives 
the continued fractions /i+ = —q/{qh'^^.^ + dn), and 
/i+ = —q/{qh^j^i + dn) with dn — (2n + v^ — a. In 
practice to find numerical values of the coefficients c„ we 
set ft,+ = and /i^j = for some, sufficiently large m 
and calculate h^ and h~ up to n — 1. Characteristic 
exponent has to determined from Eq. (jA4|) with n — 0. 

Asymptotic behavior of M^{z) and T^{z) for large z 
follow immediately from asymptotic expansion of Bessel 



functions 




M,{z) 


.^j2e-^-f^ 
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where Si, 


= Z^n=-oo''n' 



s^cos [e^^/q- 



(A5) 



,,sm{e'^-^v-l) (A6) 

,Cn{v)- To obtain asymptotic be- 
havior for large and negative z one has to join solutions 
M^(z) and T,j{z), with another pair of solution M^{~z) 
and T^(-z) at z = ^. This yields 



M,(z) 
TAz) 



z — >■ — oo 



2e^ 
7rgi/4 



mSi,cos(yge"^ + fi/- f) (A7) 
sin(^e-^ + fi.-f) 



TT (7I/4 ?7l 



cot ■nvijn — 1) cos (y^e 



7r , 



7r ^ 



(A8) 



where m = lim2^o+ My{z)/M-^{z). 



Finally we write the wave function 'l'(r) in the form 



^(r) =sin( 



R*^ L VR* 

AU In —^ 
r \ Vkr 



l''+l)\ —T.\h^ 



VR* 

\fkr 



(A9) 



where (/) is a parameter which appear in the small r ex- 
pansion Q. Now, the behavior of ^(r) at small and 
large distance described by Eqs. I11-I0J, can be readily 
obtained from asymptotic expansions HA5|I - HA8|I . 
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